Entanglement in the One-dimensional Kondo Necklace Model 
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We discuss the thermal and magnetic entanglement in the one- dimensional Kondo necklace model. 
Firstly, we show how the entanglement naturally present at zero temperature is distributed among 
pairs of spins according to the strength of the two couplings of the chain, namely, the Kondo 
exchange interaction and the hopping energy. The effect of the temperature and the presence of 
an external magnetic field is then investigated, being discussed the adjustment of these variables in 
order to control the entanglement available in the system. In particular, it is indicated the existence 
of a critical magnetic field above which the entanglement undergoes a sharp variation, leading the 
ground state to a completely unentangled phase. 
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I. INTRODUCTION 

A striking feature in quantum mechanics is the pres- 
ence of non-local correlations among different parts of 
a system. The physical property behind these purely 
quantum mechanical correlations is the so-called entan- 
glement, which has been shown to be a key ingredient 
in quantum computation and communication In- 
deed, the search of experimental proposals to implement 
a quantum computer has strongly motivated the study 
of entanglement in condensed matter systems. In partic- 
ular, schemes for quantum information processing based 
on Ising |2(, XY |3(, and Heisenberg 4] interactions have 
been suggested, leading to an intensive discussion about 
the thermal and magnetic entanglement in these mod- 
els H, IE IE @- Furthermore, the properties of en- 
tanglement have also been analysed in other condensed 
matter [ToL Hl| and interacting quantum systems 01 ■ 

Recently, initial experiments towards the construction 
of quantum gates in rare-earth-ion doped crystals have 
been carried out . The rare-earth compounds present 
several attractive properties for quantum computation 
due to the existence of a partially filled /-shell, which 
is shielded from the environment by the outer electrons. 
Such partially filled inner shell is a fundamental feature 
appearing in the heavy fermion systems |14| . character- 
ized as metallic hosts in which magnetic impurities are 
immersed. Thus, the discussion of the entanglement in 
theoretical models describing the heavy fermion physics 
turns out to be a subject which deserves some atten- 
tion. A standard hamiltonian for this purpose is the 
Kondo lattice model [lil Il7j . which assumes the pres- 
ence of one localized impurity spin on each site, cou- 
pled to the conduction electrons. From the viewpoint of 
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FIG. 1: Schematic arrangement of the KN chain with four 
sites. The spin-1/2 operators r Q and S a , where a = x,y,z, 
denote the conduction eletrons and the spins of the local mo- 
ments, respectively. 



the magnetic degrees of freedom, an important feature 
of the model is the competition between the intra-site 
Kondo screening of the localized magnetic moment and 
the inter-site Ruderman-Kittel-Kasuya-Yosida (RKKY) 
interaction among the impurity spins, mediated by the 
conduction electrons. The Kondo coupling favors a 
non-magnetic ground state, while the RKKY interaction 
tends to establish a long range magnetic order. A scheme 
for entangling two conduction electrons in the presence 
of a magnetic impurity has been proposed by considering 
the Kondo hamiltonian in ref. [lfj. 

A simplified version of the Kondo lattice hamiltonian 
was introduced by Doniach in ref. . This model, called 
Kondo necklace (KN), replaces, into the Kondo interac- 
tion, the spins of the conduction electrons by a set of 
pseudo-spins on a linear lattice, the charge degrees of 
freedom being frozen out. In spite of this approxima- 
tion, the interplay between the RKKY interaction and 
the Kondo mechanism remains as an essential feature. 
The one-dimensional KN model is defined by the hamil- 
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tonian 

N N 
H = WY^T^ + T?Tf +1 ) + J^SiA, (1) 

i=l i=l 

where rf and S" are independent spin-1/2 operators 
acting on a site i, which are given by a a /2 , with a a 
(a = x,y,z) denoting the Pauli matrices. The spin op- 
erators t° and S" are associated with the conduction 
electrons and the innercore spins hanging from the r- 
spin chain, respectively. The hamiltonian Eq. de- 
scribes a linear KN lattice with N sites, where periodic 
boundary conditions are adopted, i.e., rfj +1 = rf and 
T^ +1 = t\. The positive parameter W represents the 
hopping energy and J is the Kondo exchange coupling, 
which can be either ferromagnetic (J < 0) or antiferro- 
magnetic (J > 0). A schematic arrangement of the chain 
is displayed in Fig. Note that, since the KN model 
is defined with two spins on each site, the smallest rep- 
resentative cell of the system is described by four spins, 
each spin being taken as a quantum bit (qubit). A fully 
anisotropic version of the model, which we denote by X- 
KN, will also be considered. In this case, the band of the 
conduction electrons is represented just by an Ising term 
W J2( T i T i+i)- The isotropic KN model, given in Eq. QJ, 
will be referred from now on as XY-KN. 

The aim of this work is to discuss the thermal and 
magnetic entanglement present in the one-dimensional 
KN model. The paper is organized as follows. In section 
II, we describe the natural thermal entanglement for both 
XY-KN and X-KN models with two sites (four qubits). 
In section III, the entanglement in a larger chain is anal- 
ysed by considering the case of four sites (eight qubits). 
It is observed that, qualitatively, the entanglement of the 
system keeps the main characteristics of the case with 
only two sites. Section IV is devoted to the study of 
the effect of a magnetic field tranversely applied to the 
r-chain. Finally, in section V we summaryze our main 
results and present the conclusion. 

II. THERMAL ENTANGLEMENT IN A 
TWO-SITE KN RING 

In this section we study the entanglement between any 
two qubits of the KN model with two sites at both zero 
and finite temperature. Before describing our results we 
present the definition of concurrence, which is the mea- 
sure of entanglement used throughout this paper. The 
concurrence Q2 for a pair of qubits labelled as 1 and 2 
is defined as |lj| 

Cxi =max(A 1 - A 2 - A 3 - A 4 , 0) , (2) 

where the Xi are the square roots, in decreasing order, of 
the eigenvalues of the operator 

R = Pl 2 (a y (g) Oy)p\ 2 {(Jy <g> <Jy). (3) 



In Eq. 10, pi2 denotes the density matrix, which can 
be either pure or mixed, for the pair of qubits 1 and 2, 
and p\ 2 its complex conjugate in the the standard basis 

{| + +), | H — ), | — h), | )}. In a system with more 

than two qubits, pi 2 is obtained by tracing the density 
operator over the other qubits. The concurrence ranges 
from 0, implying an unentangled state, to 1, in which the 
two qubits are maximally entangled. 

In pure states of N qubits, it has been conjectured 
that the entanglement is distributed following the in- 
equality for the squared concurrence 

C12 + C?3 + ••• + C 1N < C^ jj, (4) 

where 6^(23. .tv) is the single-qubit concurrence, defined 
as the concurrence between the qubit 1 and the rest of 
lattice (23... N). This quantity represents the collective 
contributions of the entanglement between the qubit 1 
and all the other qubits of the system and can be obtained 
from C 1 (23..at) = 2^/det pi, where p\ is the density matrix 
for the qubit 1. A relevant aspect of the single-qubit 
concurrence is that, as it has been suggested in ref. y|, it 
may be a useful tool to identify a quantum critical point 
in a lattice system. 

A. XY-KN MODEL 

Let us consider the entanglement for the pairs of qubits 
of the two-site XY-KN model, whose spins are labelled 
with Ti, t 2 , Si, and S 2 , as represented in Fig.^ For sim- 
plicity, we will refer to T\, t 2 , Si, and S 2 , as qubits A, B, 
C, and D, respectively. The qubits A and B form a ring 
coupled via antiferromagnetic XY interaction (W > 0), 
while A and C, as well as B and D, are coupled via 
Heisenberg interaction, which can be either ferromag- 
netic or antiferromagnetic. The symmetry of the lattice 
ensures that the entanglement is the same for the follow- 
ing pairs: AC and BD, AD and BC, and AB and CD. 
Therefore we will only refer to the pairs AB, AC, and 
AD along the paper. 

The hamiltonian for the two-site model is obtained 
from Eq. Q by setting N = 2. The ground state of 
the system with either J < or J > is nondegenerate 
and can be written as the following pure state of ABCD: 

\0) xy = N xy ( I + + - -) + I - - + +) 

+ ori(| + --+> + | - + + -)) 

+ "2(| + - + -> + |- + -+»), (5) 

where the positions in the kets denote the qubits A, B, C, 
and D in this order. The normalization constant N xy and 
the functions oti and a 2 , which depend on the parameters 
J and W , are found to be 

J + 2\ X y \ X y + A^ y J — 3J^ / 4 

ai = ' " 2 = wj ' 
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y/2{l + a\+al) ' 
with \ xy being the ground state energy 



(6) 



Ky = --~2 y/Q COS - , 

. / -J(9W 2 - 16 J 2 ) ' 



54 VQ 5 



Q 



4 J 2 + iW 2 
9 



(7) 



Firstly we shall discuss the natural entanglement at tem- 
perature T=0, for which only the ground state is pop- 
ulated. If the XY coupling is taken as zero (W = 0), 
the system becomes two decoupled Heisenberg dimers, 
namely, AC and BD. In this regime, all the pairs of 
qubits are completely unentangled for J < 0. In fact, 
the absence of entanglement between qubits coupled by 
ferromagnetic Heisenberg interaction is a general prop- 
erty already pointed out in ref. 0. On the other hand, 
if J > then the qubits in the dimers are maximally 
entangled for any value of the parameter J, with the 
concurrence vanishing for all the other pairs. 

Consider now the parameter W > 0. We obtain, in the 
limit J — 0, the two-qubit antifcrromagnetic XY ring in 
the presence of two other non-interacting spins (C and 
D). In this case, the pair of qubits AB is maximally 
entangled for any W and the concurrences for the other 
pairs are zero. For the Heisenberg coupling J > 0, there 
is a competition, governed by the values of J and W, be- 
tween the entanglement of the pairs AB and AC coupled 
by Heisenberg and XY interactions, respectively. The 
concurrences Cab and Cac are shown in Figs. El and El 
the pair AD being always unentangled. By fixing the 
value of W, it is apparent that a large value of J favors 
the entanglement between AC, while small J tends to 
increase the entanglement between the qubits A and B. 
In fact, the concurrences AB and AC are complemen- 
tary, i.e., if the qubit A is partially entangled with B, 
then A can have only a limited quantity of entanglement 
with C. These results satisfy completely the inequality 
ijljl , with the upper bound limit for the entanglement be- 
ing obtained from the single-qubit concurrence between 
A and BCD, given by Ca(bcd) = 1 independently of 
the values of J and W . 

For the case J < 0, the entanglement is distributed be- 
tween the pairs AB and AD, as plotted in Figs.EJandEl 
We can see that, although there is no concurrence be- 
tween qubits coupled with Heisenberg interaction, an in- 
crease of the value of J decreases the entanglement of 
the pair AB and increases the concurrence for AD. In 
fact, it can be observed from Eqs. (|5I7|I . that a variation 
in J affects Cab and Cad, since the functions a\ and 
ct!2 depend on this parameter. As in the case J > 0, the 
inequality 10} is satisfied with C 
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= 1. 
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FIG. 2: Concurrence between the qubits A and B in the 
two-site XY-KN model for antiferromagnetic Heisenberg in- 
teraction at zero temperature. 
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In order to take into account the effect of the tem- 
perature on the entanglement in the XY-KN model, we 



FIG. 3: Concurrence between the qubits A and C in the 
two-site XY-KN model for antiferromagnetic Heisenberg in- 
teraction at zero temperature. 



consider the density matrix for the system at thermal 
equilibrium as given by p(T) — exp(~ H / kT) / Z , where Z 
is the partition function and k the Boltzmann's constant. 
In Figs. HO and we show the concurrences for the pairs of 
qubits AB and AC as functions of the temperature for an 
antiferromagnetic Heisenberg interaction and with fixed 
W = 1. By comparing the Fig. El to © and Fig. El tod 
we can see that, for the regime of very low temperatures, 
the distribution of entanglement for the pairs of qubits 
at T = is recovered. As the temperature increases, the 
entanglement for the pair AB increases, for all values of 
J, before beginning to decrease. This is explained by 
observing that, for J = 0, the ground state is degener- 
ate, being a statistical mixture of four states with energy 



FIG. 4: Concurrence between the qubits A and B in the two- 
site XY-KN model for ferromagnetic Heisenberg interaction 
at zero temperature. 




FIG. 5: Concurrence between the qubits A and D in the two- 
site XY-KN model for ferromagnetic Heisenberg interaction 
at zero temperature. 

eigenvalue — W. In this case, the entanglement AB as- 
sumes the maximum value 1 at T = 0, independently 
of the parameter W. When J is "turned on", the de- 
generacy splits up and the ground state turns out to be 
given by Eq. JSJ. However, for a nonvanishing value of 
the coupling J, there are some excited states contribut- 
ing with a greater portion of entaglement AB than the 
ground state. Thus, since the effect of the temperature 
is to populate the excited states, the entanglement AB 
turns out to increase before being completely destroyed. 
In contrast, the entanglement for the pair AC always de- 
creases. We can understand this behavior by noting that, 
for W = 0, the ground state of the system is nondegen- 
erate and the pair AC is maximally entangled for any 



FIG. 6: Thermal concurrence in the two-site XY-KN model 
with antiferromagnetic Heisenberg interaction for the pair 
AB. We have set the parameter IT to 1. 

value of J at T = 0. For W nonvanishing, the ground 
state of the system remains nondegenerate and there are 
no excited states contributing with a quantity of entan- 
glement greater than the ground state. Therefore, the 
mixing of Eq. JSJ with higher energy levels, in this case, 
can only destroy the entanglement. Moreover, as greater 
is the coupling J, greater has to be the temperature to 
lead Cac to 0. 

For ferromagnetic Heisenberg interaction, the thermal 
concurrence is different from zero for the pairs AB and 
AD. For AB it is possible to increase the entanglement 
with the temperature. However, the quantity of entan- 
glement generated is very small compared to the case of 
the antiferromagnetic Heisenberg interaction displayed in 
Fig. El 

B. X-KN MODEL 

In the anisotropic version of the KN model, the r-chain 
displayed in Fig. is composed by a set of qubits cou- 
pled via Ising interactions. As it is well known 0,0,0], 
there is no entanglement naturally present for any pairs 
of qubits in the pure Ising chain. However, when a spin- 
Si is coupled to a spin-r^ of the Ising chain, by means of 
an antiferromagnetic Heisenberg interaction, the pair of 
qubits TiSi becomes immediately entangled. It is worth 
reminding that ferromagnetic Heisenberg coupling can- 
not induce any concurrence in the system. In order to 
describe in detail the entanglement in the two-site X-KN 
model, let us consider four qubits A, B, C, and D, firstly 
at zero temperature. Similarly as defined for XY-KN 
model, the pair AB is coupled by an Ising interaction 
W > and the pairs AC and BD by an antiferromag- 
netic Heisenberg interaction J. The ground state of the 
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system can be written as 

|0)« = N x ( | + + + +) + | ) 

+ + - -> + !- - + +» 
+/?(| + - -+> + !- + + -» 

-(! + - + -) + !- + -+))) 



(8) 



where the function (3 and the normalization constant N x , 
which depend on the parameters J and W, are found to 
be 



= 



2 J + V4 J 2 + W 2 
W 



N, 



v/4(l + /? 2 



(9) 



From Eqs. (JSJ) and JJjJ, it can be shown that the entangle- 
ment for both pairs AB and AD is zero. The concurrence 
for the pair AC can be expressed as 



Cac 



B 2 



1 



(3 2 + l 



(10) 



reaching the maximum value in the limit W = 0. In 
Fig. 03 the concurrence Cac is plotted as a function of 
the parameters J and W . As it can be observed, although 
the coupling W cannot induce any entanglement for the 
Ising coupled pairs, it contributes to decrease Cac- 

The single-qubit concurrence Ca(bcd), between A and 
the rest of the chain (BCD), has also been computed. 
For J = 0, we obtain Ca(bcd) = and, for any nonvan- 
ishing J > 0, we find Ca(bcd) = 1 independently of the 
parameter W. This result indicates that, in J = 0, the 
system undergoes a quantum phase transition H char- 
acterized by a jump in the entanglement of the qubit A 
with BCD, from to 1, for an infinitesimal increase of J. 
In other words, it happens a fundamental modification in 



the pattern of entanglement present in the ground state, 
which changes from a completely unentangled state to 
one in which the upper bound limit for the concurrence 
between any two qubits assumes the maximal value. This 
quantum phase transition is very similar to that found 
in ref. 5\ for the transverse magnetic field Ising chain, 
with the coupling J playing the role of the transverse 
field. Indeed, in the presence of random coupling con- 
stants, it has been pointed out in ref. [2l| a parallelism 
of the critical magnetic properties between the X-KN and 
transverse Ising models. The results found here for the 
X-KN system indicate, in the pure case (fixed J and W), 
a similarity in the behavior of the entanglement of both 
models. 

For completeness, we show in Fig. [5] the entanglement 
of the pair AC as a funtion of temperature, where we 
have set the parameter W — 1. We observe that the 
temperature tends always to decrease the entanglement. 



III. THE ENTANGLEMENT FOR A LARGER 
CHAIN 

In this section, we numerically describe the entagle- 
ment in the KN model with a larger number of sites. By 
considering the chain with four sites, as represented in 
Fig. we show that the results found for the entangle- 
ment are qualitatively very close to that obtained earlier 
for a two-site X-KN and XY-KN ring. In fact, this seems 
to rest on the competition between the two couplings of 
the theory, which will take place independently of the 
number of sites considered. 

Here we discuss the four-site XY-KN model with an- 
tiferromagnetic Heisenberg coupling at T = 0. In this 
case, the concurrence is distributed among the Heisen- 
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1 1 

FIG. 9: Thermal concurrence for the pair of qubits AC in 
the two-site X-KN model. We have set the parameter W to 
1. 
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FIG. 10: Concurrence at T = for neighbor r-spins in the 
XY-KN model with two and four sites. We have set W = 1 
and J > 0. 



berg coupled qubits and the neighbor and next-neighbor 
pairs in the r-chain. In Figs. HOIandllllwe compare the 
entanglement by taking the chain with two and four sites. 
We call Cxy the concurrence between neighbor spins on 
the r-chain, and Ch, the concurrence for neighbor spins 
coupled by Heisenberg interaction. We observe that, for 
a fixed coupling W, the concurrence Cxy decreases faster 
with the rise of J as the number of sites becomes larger. 
On the other hand, Ch is favored if many sites are taking 
into account, increasing faster with J as we increase the 
size of the chain. The single-qubit concurrence remains 
1 and the inequality 10} is still obeyed. 

In the case of a chain with an odd number of sites 
we found that, due to the breaking of the ground state 
translational invariance [2^], the entanglement between 
neighbor qubits on the T-chain can be different from each 
other. However, this difference in the behavior of odd 
and even number of sites is expected to disapear as the 
size of the chain becomes very large, as discussed for the 
Heisenberg model in ref. 0. 
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FIG. 11: Concurrence at T = for the pairs of qubits coupled 
by antiferromagnetic Heisenberg interaction in the XY-KN 
model with two and four sites. We have set W = 1. 

IV. MAGNETIC FIELD ALONG THE Z AXIS 

We investigate here how an external magnetic field ap- 
plied in the KN model along the z direction can modify 
the entanglement of the system. The hamiltonian for this 
case is given by 

N 

H B = H + Bj2(s Z i+r*), (11) 
»=i 

where H is the hamiltonian previously defined in Eq. 
and B is a constant magnetic field. 

In order to describe the magnetic entanglement in the 
two-site XY-KN and X-KN models, we define qubits A, 
B, C and D in the same way as before. As discussed in 
section III, the qualitative results found for the two-site 
system hold for a larger chain, leading us to be concerned 
here in the former case. We first consider the magnetic 
entanglement in the KN model with antiferromagnetic 
Heisenberg coupling at zero temperature. For B = 0, 
the ground states of the XY-KN and X-KN models are 
given by Eqs. J5J and ©, respectively. The single-qubit 
concurrence in both cases is always 1, independently of 
the values of the parameters J and W. After the external 
magnetic field B is " turned on" , the ground state of the 
XY-KN model suddenly changes, as B crosses a specific 
value, from Eq. JSJ) to 

10)^ = ^(71(1 +)-|- " + -)) 

"72(1 - + --) + | + ))) (12) 

where 7V,£, is the normalization constant and the 7; are 
functions of the parameters J and W. In this state the 
single-qubit concurrence has a little decrease from 1 to 
a value which depends on J and W. For example, when 
J = W — 1, the value of the single-qubit concurrence 
is 0.989. By increasing the magnetic field to a critical 
value B c = B C (J,W) the single-qubit concurrence van- 
ishes, since the ground state becomes the fully unentan- 
gled state | ). Hence, for any fixed value of the 
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FIG. 12: The single-qubit concurrence for both two-site XY- 
KN and X-KN models with antiferromagnetic Heisenberg in- 
teraction at T=0. The parameters J and W are set to 1. For 
the XY-KN model the system presents a quantum phase tran- 
sition at B c ps 1.707. For the X-KN model the concurrence 
decays analytically with B. 
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FIG. 13: Thermal concurrence for the pairs of qubits coupled 
by antiferromagnetic Heisenberg interaction in the two-site 
XY-KN model. The parameters J and W are set to 1. As B 
crosses the critical value B c ~ 1.707 the concurrence increases 
with the temperature before vanishing at T « 0.79. 



parameter W, we obtain a critical line B c = B C (J). Tak- 
ing J = W, we can numerically interporlate the straight 
line B c s» 1.707 J. We observe that the effect of the ex- 
ternal magnetic field is the gradual alignement of all the 
spins in the chain, leading the ground state always to an 
unentangled state for a high value of B. 

The behavior for the X-KN model in the presence of 
a magnetic field is quite different to that of the XY- 
KN model. When we "turn on" the magnetic field, the 
ground state given by Eq. (jSJ generalizes to 

io)f ^i j ( ei | + + + +)+6 2 | ) 

+ e 3 | + + - -> + e 4 | - - + +> 
+ e 5 (| +__+) + | - + + _)) 

+ e 6 (| + - + -> + |- + -+»), (13) 

where is the normalization constant and the q are 
functions of the parameters J and W. By increasing B, 
all the functions e.; decrease except e 2 , which becomes 
higher. In the limit B — > oo the ground state becomes 

the unentangled state | ). In Fig.^|it is exhibited 

the single-qubit concurrence for the XY-KN and X-KN 
models at zero temperature, with the parameters J and 
W set to one. As it can be seen, the single-qubit concur- 
rence is an analytical function of the magnetic field B for 
the X-KN model, showing a smooth decrease of the en- 
tanglement of the ground state. This behavior is clearly 
different from the XY-KN case, in which the concurrence 
undergoes a sharp variation. 

Let us now consider the thermal entanglement in the 
model by taking into account the presence of the mag- 
netic field B. As discussed in Section II A for the XY-KN 
chain with antiferromagentic Heisenberg interaction, the 
increase of the temperature can rise the entanglement of 
AB due to the mixing of the ground state with some 
higher entangled excited states. In contrast, the concur- 
rence for the pair AC is always decreased as the temper- 



ature is "turned on". However, when a magnetic field B 
greater than the critical value B c is applied in the system, 
the concurrence for the pair AC can considerably be in- 
creased as the temperature rises, as exhibited in Fig. ED 
This is because, as B crosses B c , the ground state of the 
two-site XY-KN model suddenly changes from a highly 
entangled to an unentangled state. Hence, the increase 
of the temperature can here generate entanglement also 
for the pair AC due to the mixing of the unentangled 
ground state with some entangled excited states. For 
B > B c , as B approaches the critical field B c , the en- 
tanglement becomes higher with the increase of the tem- 
perature. Moreover numerical simulations point to the 
existence of a temperature for which the entanglement 
is always destroyed indenpendently of the magnetic field 
applied. For J — W = 1, this temperature is T c w 0.79 
when we compute the concurrence for a Heseinberg cou- 
pled pair. 

We have also studied the behavior of the concurrence 
in the case of ferromagnetic Heisenberg coupling in the 
presence of the transverse magnetic field. The results 
found for the single-qubit concurrence are qualitatively 
the same as the antiferromagnetic case in both X-KN and 
XY-KN models, the difference being that now a smaller 
value of B is necessary for destroying the entanglement. 

V. CONCLUSION 

In this paper we have discussed the thermal and mag- 
netic entanglement in the one-dimensional KN hamil- 
tonian, which is a simple model for describing heavy 
fcrmion systems. Concerning the natural entanglement 
present in the model, we have found that the concurrence 
for a pair of qubits is distributed according to the funda- 
mental couplings of theory. Furthermore we have been 
able to indicate, in the fully anisotropic version of the 
system, i.e. the X-KN model, a quantum phase transi- 
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tion for the Heisenberg interaction J = 0. The critical 
behavior of the entanglement has been characterized by 
means of the single-qubit concurrence, which is the upper 
bound limit for the entanglement present in the chain. 
The effect of the temperature and the presence of an ex- 
ternal magnetic field applied in the z direction has also 
been studied. These variables have been shown to be 
adjustable to control the entanglement available in the 
system. In particular, it has been found, in the XY-KN 
model, a critical magnetic field B C (J, W) which separates 
phases exhibiting a completely different pattern of entan- 
glement. 

It is known the KN model presents a rich phase dia- 
gram with remarkable critical magnetic properties [2^j |. 
It would be interesting to search a connection of this criti- 
cal behavior with the results found in this paper. Efforts 
in this direction would contribute to a more complete 
understanding of the relationship between entanglement 
and quantum phase transitions in condensed matter sys- 
tems as originally proposed in refs. 0, E| . 
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